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We suggest here a novel nano-laminate, 5[Fe]/2[WxRei_x] (x=0.6-0.8), with enhanced magnetic 
hardness in combination with a large saturation moment. The calculated magnetic anisotropy of 
this material reaches values of 5.3-7.0 MJ/m 3 , depending on alloying conditions. We also propose 
a recipe in how to identify other novel magnetic materials, such as nano-laminates and multilayers, 
with large magnetic anisotropy in combination with a high saturation moment. 



I. INTRODUCTION 

The search for new materials with tailored proper- 
ties for specific applications is a vital effort in the 
realization of novel and improved technologies. Re- 
cent developments in the ability to perform theoretical 
simulations with predictive power are just starting to 
show potential in this regard. An example here is the 
predicted huge tunneling magneto-resistance (TMR) of 
Fe/MgO/Fe trilayersf^ which subsequently was verified 
experimentallyP Most read-heads in modern disc-drives 
actually make use of this huge TMR effect. 

Magnetic materials are found in many applications. 
An example are permanent magnets, where the large en- 
ergy product of Nd2Fei4B is utilizedP e.g. in electri- 
cal machines and generators. The materials properties 
which are relevant in Nd2Fei4B is the large magnetic 
anisotropy energy (MAE) combined with a large satu- 
ration moment (M s ). Recently the risk of only a limited 
access to rare-earth elements has been discussed, with a 
potential threat for manufacturing permanent magnets 
like Nd2Fei4BP Other classes of materials, which do not 
contain rare-earth elements, but with similar magnetic 
properties, are hence in high demand. 

To identify materials with a large saturation moment 
only, is a subject of interest in its own right, with sev- 
eral challenges, e.g. as discussed recentlyP The wish of 
combining a large saturation moment with a high MAE 
becomes even more complicated, since the physical in- 
teractions responsible for one property do not necessar- 
ily correlate with those determining the other. To date 
no material has been identified with magnetic properties 
superior to rare-earth based compounds, where SmCos 
and Nd2Fe 14 B stand out in particular. Several new ma- 
terials have however been considered, where the tetrag- 
onal compound FePt is one of the most recently studied 
materials P 

Theory based on density functional theory has proven 
accurate in reproducing measured saturation moments, 
and the difference between observations and measure- 
ments is seldom larger than 5-10 %. As regards the 
magnetic anisotropy it is significantly more difficult, 
since the accuracy needed is very high. However, some 



developments have recently been madeP As an exam- 
ple we mention the predicted large MAE of a tetrago- 
nally strained FeCo-allojP which subsequently was veri- 
fied experimemtallyP^l 

In this communication we use first principles theory 
and focus on the magnetic anisotropy energy (MAE) and 
suggest a novel nano-laminate with enhanced magnetic 
hardness in combination with a large saturation moment. 
In addition to identifying this new material, we propose 
here a recipe in how to identify other novel magnetic 
materials, such as nano-laminates and multilayers, with 
large magnetic anisotropy in combination with a high 
saturation moment. 



II. GEOMETRY AND COMPUTATIONAL 
DETAILS 

In a similar spirit as that used in Ref@]we wish to inves- 
tigate the MAE as a function of electron concentration, 
but in the present study it is the electron concentration 
of the ligand atoms, W and Re, which are under focus, 
not the atoms carrying the magnetic moment (Fe). The 
advantage of tuning the electron concentration of the lig- 
and states, is that they can be chosen with a large spin- 
orbit coupling, which via hybridization with the Fe 3d 
states are expected^ to influence the anisotropy of the 
whole multilayer. We considered multilayers of differ- 
ent thickness of the Fe and W-Re layer, and the chosen 
composition, 5[Fe]/2[W x Rei_ I ] (x=0.6-0.8), is the one 
that gives a suitable strain of the W-Re layer, as well as 
an optimized electron concentration, as described below. 
The final geometry was obtained via force minimization 
using the VASP package. This resulted in essentially a 
tetragonal structure with an effective c/a ratio of the W- 
Re layers (and of the Fe layers), i.e. the ratio between 
atomic distances of the W-Re layer in the out-of-plane 
and in-plane direction. As we shall see below this ratio 
is a crucial parameter determining the MAE. 

Using the geometry optimized structure we then used 
two all electron methods to calculate the MAE, which are 
described in detail below. The two methods are below re- 
ferred to as FPLMTO-VCA^and KKR-CPAPln these 
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calculations we also varied, within limits, the c/a ratio 
of the W-Re layer, and as we shall see below, this causes 
drastic changes in the MAE. An experimental realiza- 
tion of modifying the effective c/a ratio of multilayers is, 
in general, readily done by varying the thickness of the 
individual layers. 

As mentioned, from the geometry optimized structure 
we used an all electron, full-potential linear muffin-tin 
orbital (FPLMTO) method^ to evaluate the MAE. This 
calculation used the simple virtual crystal approximation 
(VCA) to treat the alloying between W and Re. In these 
subsequent calculations we also varied, within limits, the 
c/a ratio of the W-Re layer, and as we shall see below, 
this causes drastic changes in the MAE. An experimental 
realization of modifying the effective c/a ratio of the W- 
Re layers can be done by tuning the relative amount of 
Fe and W-Re layers, i.e. by considering the superlattice 
m[Fe]/n[W x Rei_a;], with m/5 and n / 2. 

In order to give further credit to our theoreti- 
cal study we repeated the MAE calculations for the 
5[Fe]/2[W a; Rei_ a ;] superlattice in terms of the fully 
relativistic screened Korringa-Kohn-Rostoker (SKKR) 
method 13 . This method can combine a full relativistic 
description with the more advanced coherent-potential 
approximation (CPA) to tackle the chemical disorder in 
the W x Rei_ x layers (where x=0.79). We used the lo- 
cal spin-density approximation of the density functional 
theory as parametrized by Vosko et aZPSI, the effective po- 
tentials and fields were treated within the atomic sphere 
approximation (ASA) with an angular momentum cut- 
off of i m ax = 3. To calculate the MAE we used the 
so-called magnetic force theorem^ in which the previ- 
ously determined self-consistent effective potentials and 
fields are kept fixed and the change of total energy of 
the system with respect to the direction of the magne- 
tization, n, is approached by that of the single-particle 
(band) energy. In these calculations the energy integra- 
tions were performed by sampling 20 points on a semi- 
circular path in the upper complex semi-plane and, for 
the energy point closest to the real axis, nearly 10.000 
/c-points were selected in the surface Brillouin zone that 
ensured a numerical error not exceeding about 5% of the 
calculated MAE values. 



III. RESULTS 

We show in Figjlji the calculated MAE of 
5[Fe]/2[W a; Rei_ 2 ;] as a function of the local c/a ra- 
tio of the W-Re layers. The starting concentration 
which was considered was x=0.8. In the figure we show 
both results from FPLMTO- VCA as well as KKR-CPA. 
From the figure it is clear that the two methods give 
similar overall behavior, but that the KKR-CPA calcu- 
lations result in a maximum of ~ 4 meV/f.u., whereas 
the FPLMTO- VCA results have a maximum of - 3 
meV/f.u. This is comforting since the KKR-CPA treats 
the disorder better while the FPLMTO- VCA handles 



the relaxed structure more accurately. The two methods 
give both a similar shape of the MAE curve with respect 
to c/a ratio of the W-Re layer, and the maximum of 
the curves occur at similar c/a. The maximum value 
of the KKR-CPA results can be compared to the MAE 
of FePt, which from theory is ^3 meV/f.u.p^ whereas 
experiment results in 1.2 meV/f.uP^ The MAE of 
5[Fe]/2[Wo.8Reo.2] is hence very comparable to that of 
FePt. A further tuning of the magnetic anisotropy can 
be done by modifying the amount of alloying between 
W and Re in the multilayer, while keeping the c/a ratio 
the same as that giving the maximum MAE value in 
Figjl] i.e. 1.34. Calculations based on FPLMTO- VCA 
are shown in FigJlJ), and we note that the MAE can be 
made ~ 30 % larger by tuning x to 0.6. 

As to the magnetic moments they are rather constant 
for the different values of x considered, and we obtain 
values of 1.65 \ib per interface atom and 2.4 fig for the 
sub-interface Fe atoms, and 2.2 fis for the Fe atom in 
the middle of the Fe slab. The moment on W-Re is -0.1 
fiBi resulting in a total moment just above 10 /xg/f.u. 
For comparison we note that for FePt this value is 3.3 
Mb/Lu\ 

So far we have only considered calculated values 
per formula unit, whereas for practical applications 
it is with respect to volume one should compare the 
5[Fe]/2[W a; Rei_ a ;] system with e.g. FePt. As concerns 
the MAE we obtain a maximum value 5.3-7.0 MJ/m 3 
(depending on concentration or computational method) 
which should be compared to a measured value of ~ 
7 MJ/m 3 for FePiR For the moments we obtain 0.11 
^«b/A 3 which is to be compared to 0.12 /x^/A 3 for FePt. 
A final note on the comparison to the properties of FePt, 
is that this material often exhibits partial chemical dis- 
order on the Fe and Pt sublattices, which drastically re- 
duces the MAEpl 

In order to obtain a more detailed, microscopical un- 
derstanding for why the MAE of Fig{l] obtains a maxi- 
mum at c/a close to 1.3 we show first of all that the large 
MAE is caused by the large spin-orbit coupling of the 
W-Re atoms. As a second step we will illustrate that it 
is a delicate detail of the energy bands of the W-Re layer, 
which take a special structure for a c/a ratio close to 1.3, 
which determine the large MAE. We start however, with 
analyzing the spin-orbit effect of the W-Re atoms. 

As a direct way to illustrate the dominant role of the 
W x Rei- x layers in the formation of the high perpendic- 
ular MAE of this system, we varied the spin-orbit cou- 
pling strength at these layers for the optimal c/a ratio 
(c/a — 1.3) and concentration (for x = 0.8, other val- 
ues of x give a similar picture). Here we employed the 
technique developed by Ebert et al. within the KKR for- 
malism.^ The corresponding MAE values are shown in 
Fig. [2] for SOC strengths (£) scaled from to 1. Remark- 
ably, for £ = 0, i.e. when SOC is considered just in the 
Fe layers, the MAE is slightly negative preferring thus 
an in-plane orientation of the magnetization. Increasing 
£ in the W^Rei-^ layers the MAE starts to increase as 
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FIG. 1: (Color online) (a) Calculated magnetic anisotropy of 
5[Fe]/2[W a: Rei_ a: ] (x=0.8) as a function of effective c/a ratio 
of the W-Re layers, (b) Calculated magnetic anisotropy of 
5[Fe]/2[W a; Rei_ a; ] as a function of W concentration, x, for a 
c/a ratio of 1.34 of the W-Re layers. The MAE is shown per 
formula unit. 



a quadratic function: this behavior is well understood in 
terms of the second order perturbation theory, discussed 
below. For about 0.5 < £ < 1, however, the MAE seems 
to follow a linear behavior implying that in this regime 
SOC can not be regarded as a small second order per- 
turbation to the scalar-relativistic spin-polarized band 
structure. Overall, the data of Fig[2] demonstrate that 
the large MAE is indeed caused by the large spin-orbit 
coupling of the W-Re atoms. 
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FIG. 2: Calculated MAE (circles) within the KKR-CPA ap- 
proach with SOC strength scaled between and 1. The 
optimal setup, c/a — 1.3 and x = 0.8, is fixed for the 
5[Fe] 2[W a; Rei_ a ;] superlattice. Solid line serves as a guide 
for the eye. 



We now turn to investigating why an effective c/a ratio 
of the W-Re layer close to 1.3 causes a maximum value 
of the MAE. The MAE of a nano-composite material can 
be expressed as 

AE so = J2 AE qq , = [Etf fa) - E qq , (n 2 )] , (1) 

qq' qq' 

where q is the atomic species and hi are two spin- 
quantization axis, one out of plane and one in plane. 
Each of the E qq i (n) can now, according to second order 
perturbation theory, be written as 

E q q'{ri) ^ ^ ^ ^ ^ ^ ^ , kis,gm,q'm / ^'k^s / ,g / m",gm" 

kij {s} {m} 

(qms \H SO (n) | qm"'s'} {q'm"s'\H S o (n) \ q'm's) 
x ) 

where k represents the sampling points in the Bril- 
louin zone, i, j is the occupied and unoccupied states, 
{s} = {s, s'} run over the spin components and {to} = 
{to, to', to", to'"} run over the magnetic quantum num- 
bers. The basis functions \qlms) are specified through the 
atomic site q, the I, to and s quantum numbers, i.e. the 
azimuthal, magnetic and spin component quantum num- 
bers. Note that the band character nki S ,qm,g'm' includes 
hybridization, i.e. mixing of different basis functions at 
e.g. different atomic sites q. The contributions from Eqn. 
[2] are seen to become large if the energy difference be- 
tween occupied and unoccupied electron states is small 
to permit a strong coupling of the spin-orbit interaction, 
especially if this happens for a large fraction of the BZ. 
This fact is further enhanced if the states that couple 
via the spin-orbit coupling have a significant contribu- 
tion from the states residing on the heavy W-Re atom, 
since this guarantees that at least one of the spin-orbit 
matrix elements in Eqn. [2] is then large. 

It is in such an analysis that the effective c/a ratio of 
W-Re close to 1.3 becomes clear. Inspection of a W-Re 
alloy reveals that as a function of c/a ratio, conspicu- 
ous features of the band structure explain the enhanced 
MAE. To illustrate this we show in Fig. [3] the energy 
bands along r - M - G - P - Y - M - T - X (of the bet 
BZ) of a W-Re alloy (calculated in this case by virtual 
crystal approximation - VCA - for simplicity, and with 
x=0.8). From the figure one may note a spin-degenerate 
band situated very close to the Fermi level (E-p) between 
Mo and T and X. As a function of changing the c/a ratio 
this band moves from being completely occupied to be- 
come unoccupied, and for c/a=1.3 it is situated right at 
Ei?. When put in contact to Fe, in the 5[Fe]/2[W a; Rei_ a ;] 
multilayer, the narrow band discussed in Figj3] will be- 
come exchange split, due to the induced moment of the 
W-Re layers. In order to illustrate this effect in a simple 
way, we have performed a fixed spin moment calculation 
of W-Re, with the same induced moment on the W-Re 
site as for the 5[Fe]/2[W 2; Rei_ 2 ;] multilayer. The result- 
ing dispersion relationship is also shown in Fig|3] (bottom 
panel). Here we note a breaking of the spin-degeneracy 
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close to Ep , and the flat band is split into up- and down- 
bands, which for a substantial part of the BZ are situ- 
ated on each side of Ep in a large part of the in-plane 
Brillouin zone. Spin-orbit coupling via the Z+s_ + Ls + 
terms of H$o m Eqn.2 then provides a large contribu- 
tion to the magnetic anisotropy. A modification of the 
c/a ratio away from this value moves the flat bands away 
from Ep, which reduces the MAE. 




FIG. 3: (Color online) Non spin-polarized bands (upper 
panel) and spin-split bands (lower panel) of W-Re for bet 
c/a ratio 1.3. 



IV. CONCLUSION 

We have here presented theoretical evidence for that 
a non-laminate involving W-Re layers in contact to Fe, 
can result in a magnetic material with a substantial mag- 
netic moment and a large magnetic anisotropy. The new 
material compares in performance to the well known ma- 
terial FePt, both in terms of saturation moment as well 
as anisotropy. One advantage with the present system is 
the possibility to tune the MAE by tuning the local c/a 
ratio of the W-Re layer, which can be done by introduc- 
ing more Fe layers. This simultaneously also increases 
the saturation moment per volume. 

Our analysis points to a general route to find materials 
with a large magnetic anisotropy. A strong ferromagnetic 
material, preferably based on Fe or an FeCo alloy, should 
be in contact with a heavier material with large spin-orbit 
splitting. Most crucially, the electronic structure of the 
heavy material should be similar to the schematic energy 
band plot of Fig|4j Here bands which are dominated by 
states on the heavy atom are drawn to be exchange split, 
due to the interaction (hybridization) of the ferromag- 
netic material. If a large fraction of the BZ has majority 



bands below the Fermi level, and a large fraction of the 
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FIG. 4: (Color online) Schematic picture of the electronic 
structure of bands projected on heavy atoms, having possi- 
bility of inducing a large magnetic anisotropy. Majority spin 
states in red, minority states in blue, and the Fermi level a 
dashed horizontal line. The exchange splitting is marked as 
A*. 



BZ has minority states above (as Figj4]is drawn), then 
the spin-orbit coupling of the heavy atom will give rise to 
a large contribution to the MAE for many wave vectors 
(as Eqn.2 shows). 

The situation shown in Figj4] corresponds to a set of 
narrow bands close to the Fermi level, which corresponds 
to a system with a large value of the density of states 
at the Fermi level (DOS(Ep)). Hence, a large value of 
DOS(Ep) is also a requirement when searching for novel 
materials with a large MAE. Unfortunately a large value 
of DOS(-Ef) is often associated with an unstable or meta- 
stable crystal structure,^ which points to that it is meta- 
stable systems, such as multilayers and nano-laminates, 
which one should consider, when searching for novel ma- 
terials with a large MAE. 
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